The role of the lattice structure in determining the magnon-mediated interactions 
between charge carriers doped into a magnetically ordered background 
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We use two recently proposed methods to calculate exactly the spectrum of two spin-| charge 
carriers moving in a ferromagnetic background, at zero temperature, for three types of models. By 
comparing the low-energy states in both the one-carrier and the two-carrier sectors, we analyze 
whether complex models with multiple sublattices can be accurately described by simpler Hamil- 
tonians, such as one-band models. We find that while this is possible in the one-particle sector, 
the magnon-mediated interactions which are key to properly describe the two-carrier states of the 
complex model are not reproduced by the simpler models. We argue that this is true not just for 
ferromagnetic, but also for antiferromagnetic backgrounds. Our results question the ability of simple 
one-band models to accurately describe the low-energy physics of cuprate layers. 

PACS numbers: 71.10.Fd, 71.27.+a, 75.50.Dd 
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I. INTRODUCTION 

Many materials under current study, such as cuprates, 
manganites, pnictidcs, irridates, etc., have complicated 
structures with several types of atoms in the basis. In- 
cluding orbitals for all these atoms in a model Hamilto- 
nian would make it impossibly difficult to solve, besides 
introducing unreasonably many parameters. 

Thus, the first major challenge in studying such mate- 
rials is to understand what is the minimal model Hamil- 
tonian that properly captures their low-energy proper- 
ties. (The second major challenge, of course, is to figure 
how to solve it). Some steps in this process are fairly 
straightforward. For instance, many of these materials 
have layered structures, and there are many indications 
that the interesting physics is hosted by certain layers 
which are common to all compounds in the family. It is 
therefore an easy decision to start by modeling one such 
layer - for example, a Cu02 layer for cuprates. 

Even though this is a significant simplification, the 
proper minimal model to describe the low-energy prop- 
erties of one layer is still not obvious. To continue with 
the CuC>2 example, one could write a model that includes 
the 3d x 2_ y 2 orbitals for Cu, which are known to be half- 
filled in the undopcd parent compound, plus the appro- 
priate 2p x / y O orbitals. These are filled in the undoped 
parent compounds, but the cuprates are charge transfer 
insulators 1 and these are the orbitals expected to host 
the holes introduced by doping. This approach leads to 
the three-band model of Emery, 1 although for later con- 
venience we prefer to label it here as a two-sublatticc 
model: one for the Cu sites, effectively hosting spins- \ if 
double occupancy is forbidden, and one for the O sites, 
hosting the holes (the charge carriers). 

A simpler option is a one-band Hubbard model which 
describes effective states located on one lattice (at the Cu 
sites, for Cu02 layers), or its even simpler counterpart, 



the t- J model, which additionally discards all doubly oc- 
cupied states and describes carriers in a spin background. 
Even these simplest possible interacting Hamiltonians, 
which depend on only one dimensionless parameter (U/t 
or J/t, respectively) are still far from completely under- 
stood in the strongly correlated limit, despite a tremen- 
dous amount of effort.'' 

It is the difficulty to solve such correlated Hamiltonians 
that explains the drive to simplify the models that we 
study to the utmost possible. However, it is reasonable to 
expect that as more and more details are discarded, there 
is a point beyond which the model no longer contains 
the physics that we aim to study. Arc the Hubbard or 
t-J models already past this point, as far as cuprates are 
concerned, or do they still contain the low-energy physics 
of the more complex Emery model? 

An affirmative answer to the latter question was given 
by Zhang and Rice, 1 who also clarified the nature of the 
states appearing in these one-band models, namely the 
Zhang-Rice singlets (ZRS). We review their arguments in 
some detail below; for now it suffices to say that the ZRS 
is a composite object involving a doping hole occupying 
a certain linear combination of O orbitals, and which is 
locked in a singlet with a hole (spin) at a Cu site. 

Zhang and Rice's work on the equivalence between the 
low-energy properties of the Emery and t-J models is 
based on arguments limited to the one-carrier (one dop- 
ing hole) sector. Its validity is still not fully settled be- 
cause of the lack of exact results for these Hamiltoni- 
ans. Accurate analytical approximations are not avail- 
able, while conclusions based on numerical work are ham- 
pered by the restriction to rather small clusters. Recent 
results for clusters with 32 Cu02 unit cells suggest that 
at least in parts of the Brillouin zone (BZ), the quantum 
numbers and symmetries of the lowest eigenstate of the 
Emery model are different from those of its t-J counter- 
part, although the energy dispersion is rather similar. ' 

While efforts to understand if the ZRS offers a good 
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description of the quasiparticle are on-going, it is impor- 
tant to note that even if this is proven true, it would still 
not settle the question whether one-band models give a 
proper description of the low-energy properties of weakly 
doped CuC>2 layers. This is because a model Hamilto- 
nian must describe properly not only individual quasi- 
particles, but also their interactions. This is especially 
important in cuprates where high-temperature supercon- 
ductivity appears upon doping. Much of the commu- 
nity believes the pairing to be mediated through magnon 
exchange, 1 ' although phonon-mediatcd pairing' and even 
combined mechanisms^ are also favored. Clearly, the 
model Hamiltonian must describe properly the magnon- 
mediated interactions before one can decide if they are 
sufficient to facilitate pairing and superconductivity at 
these rather high temperatures, or whether coupling to 
additional bosons, e.g. phonons, needs to be considered. 

To settle this issue, one must compare the low-energy 
properties of the two models in the two-carrier sector; 
the single-carrier sector has no information about inter- 
actions. For the Emery and t-J models, such a compari- 
son is even more difficult, since finite size effects become 
more important as the number of holes increases. Recent 
efforts in this direction were rather inconclusive. 9 

Although the discussion so far was focused on cuprates, 
similar questions appear for other materials. Is it reason- 
able to ignore 10 the ligand O and use models for man- 
ganites which only involve Mn 3d orbitals? Do the As 
ions play any role 11 in the low-energy physics of pnic- 
tides? These and many other similar questions have, as 
a common denominator, the issue of whether one has to 
use complex models with different ions on different sub- 
lattices to properly describe such materials, or whether 
it suffices to use one-lattice models, possibly based on 
some composite states which effectively account for the 
role played by other ions, such as the ligands. 

Here we focus on one class of such questions, for ma- 
terials where the parent compound has long range mag- 
netic order and is a charge-transfer insulator. (A short 
version of this work was published in Ref. 12). Our con- 
clusions are based on results for the one- and two-carrier 
sectors for a background that has ferromagnetic (FM) 
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FIG. 1: (Color online) Sketch of a CuCh plane. If double 
occupancy is forbidden, each Cu site has a hole in a ?>d x i_ y i 
orbital, behaving like a spin-| (red arrows). Upon doping, the 
O 2p x / y orbitals (empty blue dots) accommodate the doping 
holes (filled blue dots, with arrows indicating the spins). 



order. This is much simpler than an antiferromagnctic 
(AFM) background and allows us to find exact solutions, 
so that we are certain that differences between results 
are due to the models themselves. We find that the sim- 
pler one-lattice models do not properly describe magnon- 
mediatcd interactions between carriers, and identify the 
simple reasons for this. Moreover, we argue that these 
conclusions are relevant for AFM backgrounds as well. 
We believe that our results offer strong arguments that 
the one-band Hubbard and t—J models are inappropriate 
to describe cuprate layers. 

The work is organized as follows. In Section II we 
briefly review the ZRS arguments. In Section III we in- 
troduce our models, which mirror the steps in the ZRS 
derivation, but for a FM background. Section IV sum- 
marizes the formalism (details are presented in various 
Appendixes). The results are in Section V, and Section 
VI contains a detailed discussion and our conclusions. 



II. BRIEF REVIEW OF THE ZRS 

The Emery model for a Cu02 plane, sketched in Fig. 
1, includes explicitly both the Cu and the O sublatticcs, 
with their 2>d x 2_ y 2 and 2p x / y orbitals, respectively. The 
goal is to try to reduce it to a one-lattice model, built on 
effective states located at the Cu sites. 

Since the doping hole occupies O 2p orbitals, the first 
step to remove the O sublattice is to make linear combi- 
nations of the four O orbitals neighboring one Cu, which 
are thus centered at the Cu sites i. These are described 
by the f/f! operators in Eq. (5) of Ref. 2. Because 
of phase coherence, the "S" combination where the O 
orbitals have x 2 — y 2 symmetry is the low-energy state, 

and from now on we restrict the discussion to it. 

(s) 

The P> a J states are not orthogonal since neighbor 
Cu share an O ligand, [Eq. (6) of Ref. 2]. To fix 
this, Zhang and Rice first define their Fourier trans- 
forms Pk,o- = -7^ Ac J2i e ~ lkRi Pia\ where the normal- 
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ization factor /?k = 



■i (cos fee + cosfcj,) is due to 



the non-or thogonality of the P> J [Eqs. (8) and (9) of 
Ref. 2]. It is problematic that /3k diverges near the cen- 
ter of the BZ, but ignoring this fact one can Fourier 
transform back to real space to obtain the operators 



/Ns 



ikR, 



-Pker- These are now properly nor- 



malized and describe "effective Cu orbitals", although 
they are in fact rather complicated linear combinations 
of O orbitals. The next step is to note that exchange 
favors locking a doping hole occupying such an orbital 
in a singlet with the hole (spin) located on the central 
Cu. This results in the ZRS described by the operators 
V>i = 73 (fafdu. ~ 0i|dif), see Eq. (10) of Ref. 2. Zhang 
and Rice then argue that if double occupancy is forbid- 
den, a t-J Hamiltonian based on these ZRS captures the 
low-energy states of the Emery model. 
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FIG. 2: (Clor online) Sketches of the three models. Models 
I and model II have two bands. One (empty circles) hosts 
charge carriers introduced by doping (filled circle with arrow 
for the carrier spin) and the other one hosts lattice spins (red 
arrows). In model I the carriers propagate on a different sub- 
lattice than the one hosting the spins, while in model II, they 
are on the same lattice. Model III has only one band which 
hosts both spins (arrows) and ZRS-like objects (filled circle). 



III. MODELS 

In an effort to mirror these steps in going from the 
Emery to the t-J model, we study three cases sketched 
in Fig. 2. For the sake of simplicity we restrict ourselves 
to one dimension (ID). Generalizations to higher dimen- 
sions are straightforward, but as explained below, they 
do not lead to qualitative changes. 

Case I is the "parent" two-sublattice, two-band model. 
One sublattice hosts the lattice spins, of quantum num- 
ber S, while the other hosts the charge carriers. This is 
the ID counterpart of the Emery model for the 2D Cu02 
layers if double occupancy is forbidden (apart from the 
FM order, of course). Model II is a one- lattice, two- band 
model. One should roughly think of these on-site carrier 
orbitals as the analogs of the 4>i,a m the ZRS derivation. 
Finally, model III is an even simpler one-band model, the 
equivalent of the t-J model obtained after locking doping 
holes in singlets with lattice spins. These parallels will 
be made more clear in the Results section. 

In all these models lattice spins interact ferromagneti- 
cally with their nearest neighbors: 



(i) 



where Si is the operator for the spin at site Ri = ia, and 
a = 1 is the lattice constant. Its ground-state is |FM) = 
| + S, +S, ■ ■ ■ + S) and has zero energy. For reasons de- 
tailed below, the excited states of interest will be the one- 



magnon states S q |FM) = 



energy D, q = 4J5sin 2 (|J. The raising and lowering op- 
erators Sf 1 = Sf±iSf correspond to creation (— ) and an- 
nihilation (+) operators for magnons, respectively. If we 



rewrite Hs = 
S 2 ) and H x s ' v 



H s + Hg , 

= -4v. 

2 t-Ji 



where 



that He" describes magnon hopping. 



it is clear 



-Q+ 
D i+1 



In models I and II, the second band hosts spin-i charge 
carriers (electrons or holes) introduced by doping. For 
model II these carriers propagate on the same lattice as 
the spins. For model I, they move on a different sublattice 
which is interlaced with the spin sublattice (see Fig. 2). 
For both models the carriers arc described by a Hubbard 
Hamiltonian He = T + U , where: 

T=-tY.( c l + s+i,^+s,*+n.c.), (2) 

E,tr 

U = U ^ n i+s,t n i+s,i, (3) 



for model II. Here, cL_ s 



and 6 = i for model I, 5 

is the creation operator for a carrier with spin a at site 
i + 5, and n i+ s,a = c\ +s a c i+&,cr is the occupation number 
operator. The hopping T is diagonalizcd by the states 



J 

-fe.cr 



1 V 



jk-n 



+ s c i j_ x , where TV is the number of 

2 + 0,(7 1 



unit cells. The resulting free carrier dispersion is e(k) = 
— 2£cos(fc), where k E (— n, ir] is in the BZ. 

The exchange between the carriers' spins and the lat- 
tice spins takes the simplest form: 



H 



(/) 



for model I, while for model II: 



Here, s l+s = Y^ a 



(3 C i+8,a 2 



" c i+<5,/3 is the spin operator 



for carriers at site i + 5, and a are Pauli matrices. In the 
following we choose Jo > to be AFM, as appropriate 
for a Cu02-like situation. This also has the advantage 
that it favors an infinitely lived low-energy quasiparti- 
cle, as detailed below. However, our solution works for 



FM exchange, Jq < 0, as well. If we write T-L 



(iji) _ 

ex — 



ex\z 

(in 

ex;x,y 



with, e.g., % 



(ii) 

ex\z 



4 L E t [4 s r + S 7S+], etc., the latter term 



describes spin-flip processes in which a charge carrier ei- 
ther absorbs or emits a magnon. 

While in this work we take the total Hamiltonians for 
cases I and II to be Hiji = Hc+'Hs+'Hi I x II \ it is impor- 
tant to note that if the no double-occupancy restriction 
is imposed properly, there are additional terms arising 
in the two-sublattice analog of the Emery model. 5 ' 9 ' 13 
For example, exchange between lattice spins is blocked 
on bonds where carriers reside. Since this is of order 
J, which will be chosen as the smallest energy scale 
in the problem, its effects may not be that important. 
More substantial is the appearance of a "spin-swap" 
term, which describes carrier hopping while its spin is 
exchanged with that of the lattice spin past which it 
hops. This term has an energy scale comparable to Jq, 
as expected since they both come from second-order pro- 
cesses in the t p d hopping. 9,1,5 Since Jq is a large energy 
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scale, the effect of ignoring this term may be substantial. 
Such terms can be easily handled by our exact solution, 
but we do not include them here so as to keep models I 
and II as similar as possible. This allows us to conclude 
that qualitative differences between their results are due 
to the different lattice structures, whereas if such addi- 
tional terms were present in one but not the other model, 
this would no longer be clear. Note that these additional 
terms were ignored by Zhang and Rice as well. 

The one-band model III is described by a t-J model, 
where the spin-spin exchange Us of Eq. (1) is supple- 
mented by "hole" hopping (where the "hole" is a charged 
singlet, not a regular fermionic charge carrier). Because 
all lattice spins are up, this is completely equivalent with 
a Hubbard Hamiltonian; double occupancy is automati- 
cally prevented here by Pauli's principle. 

Models I and II have three energy scales: t, J and 
Jo. We will assume throughout that J <C i, Jo- This 
is physically reasonable since this super-exchange arises 
from 4th order processes in the t p d hopping, whereas Jo 
and t (which is similar to the spin-swap term but without 
exchange of the spins) are second order in t p d- Of course, 
t could also include a direct t pp contribution. Which of 
the latter two is larger depends on the interplay between 
the Hubbard repulsion on each of the two sublattices (Ud 
and Up, in cuprate language) and the charge transfer en- 
ergy (A). In model III, the Jo energy scale has been 
integrated out. It stabilizes the low-energy nature of the 
ZRS-like object but it does not influence its dynamics, 
hence only the t and J energy scales remain. 

Finally, note that all these Hamiltonians are invariant 
to translations and commute with the z-component of 
the total spin S£ ot = + s i+s)- This has impor- 

tant consequences for the structure of their eigenstates, 
as described next. 



IV. FORMALISM 

Before we can understand what interactions arise be- 
tween quasiparticles, we need to first understand the 
quasiparticles themselves. This is achieved in the sin- 
gle carrier sector. We review the solution for this first, 
and then move on to the two-carrier sector. 



A. Single charge carrier sector 

Since the FM background breaks spin rotational sym- 
metry, we need to treat separately the cases with the 
charge carrier injected with spin up and spin down. 

The former case is solved trivially at T = 0. Due 



to the restriction to the S? 



NS 



spin- flips are possible. The state c^|FM) is therefore 
an eigenstate of with energy E k _^ = e(k) + jJqS, 

with 7 = 1 for model I and 7 = \ for model II. Model 
III cannot differentiate between an undoped system and 



a system doped with spin-up carriers, and is therefore 
unable to describe the physics of this case. 

The T = solution for a single spin-down charge car- 
rier in model II has been known for a long time, and 
has recently been generalized to two-sublattice models 
like model I. 14 An exact solution can be obtained for the 
Green's function: 



(FM|c fc4 GM4jFM), 



(4) 



where G(w) = [lo — H + 277] _1 is the resolvent for T-L, 
h = 1 and i] is an infinitcsimally small positive number 
which enforces the retardation condition. Physically, it 
corresponds to the introduction of a finite lifetime ~ 1/?/. 
Using a Lchmann representation: 



G ; (fc, W ) = ]r 



(FM|c fe4 |^)(^l4 4 |FM) 



IT) 



(5) 



we find that Gi(k,uj) has poles at the eigenenergies e n 
in the one-carrier, S£ ot = NS — \ subspace: H\ip n ) = 
EnlVVi)- The weights at these poles are the overlaps be- 
tween eigenstates \tp n ) and free-particle states c\ , |FM). 

The way to calculate Gi(k,u>) is well established. 14 ' 15 
One uses the Dyson identity: 



G(w) = G (w) + G(u)VGo(w), 



(6) 



where Go(w) is the resolvent for Ho = T + H s and 
V = % — Jio- The spin- flip part of the carrier-spin 
exchange, Wex*x}y, mixes the state c£ JFM), of energy 
Ef~ ± = e(k) — "/JqS, with the continuum of one-magnon 
states c' k ^S~\FM), of energy Ek- q .\ + £l q - Due to the 
restriction to the NS— ^ subspace, the excitation of more 
than one magnon is not allowed. The resulting two cou- 
pled equations of motion (EOM) can be easily solved 
analytically (for more details see Refs. 14,15). 

For model III, the solution in this subspace is trivial: 
a "hole" propagates freely with energy e(k). 



B. Two charge carriers sector 

Here, there are three distinct cases: (i) both carriers 
are injected with spin up, S^ ot = NS + 1; (ii) carriers are 
injected with opposite spins, S* ot = NS; and (hi) both 
carriers are injected with spin-down, S* ot = NS — 1. 

Case (i) is trivial at T = 0, since no spin-flip processes 
are possible, and therefore there is no magnon-mediatcd 
interaction. (The Hubbard repulsion has no effect, either, 
since both carriers have spin- up). The eigenstates are 



subspace, no simply c\ ^c\, JFM) with energy E kA + E k > it . 



The T = solution for case (ii) is discussed in detail 
next. We present two methods to obtain this solution, 
one based on a /c-space formulation and one based on 
a real-space formulation. The former is exemplified for 
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model II while the latter is exemplified for model I, how- 
ever both models can and have been solved with both 
methods. For model III we run into the same problem 
as before: since this model cannot distinguish between 
lattice-spins and doped spin-up carriers, we cannot con- 
sider this case for it. 

Case (iii) turns out to be less interesting for our pur- 
poses. It is briefly discussed at the end of Section V. 

1. k-space solution for model II in the two-carrier, 
Stot = NS sector 

This solution is rather similar to that for G^(k, uj). We 
start by defining the following Green's functions: 

G(k,q,q',u } ) = (k,q'\G(u>)\k,q), (7) 

where \k,q) — c' k + ^ _^ JFM) is a two-carrier state 

with total momentum k and S* ot = NS. While the total 
momentum k is conserved, q is subject to change due to 
on-site Hubbard scattering and to magnon-mcdiated in- 
teractions. The latter describe the emission of a magnon 
by the spin-down carrier and its subsequent absorption 
by the spin- up carrier. As a result, the two carriers ex- 
change their spins, besides exchanging some momentum. 



Note that we obtain only two coupled EOM, since the 
initial states are linked only to states with two spin- 
up carriers and one magnon. Furthermore, we do not 
need to know the full F to solve for G, only its av- 
erage jj J2 P F{k 7 q,q' ,p 7 uj). We therefore need to ex- 
press this average in terms of G and Go- This al- 



In order to construct the EOM for G(k,q,q' ,uj) we 
split the Hamiltonian into two parts, rl = rlo + V where 
rlo = T + ri z s , and repeatedly use Dyson's identity, Eq. 
(6). After using it once we obtain: 



G(k, q, q',uj) = y q . q , + - £ G(k, Q, q', u) 

+^E F ( fc ^,<z',P,^) 1 Go(*+**-9,w). (8) 
p ) 



Here Go(k, k',uj) = [uj — e(k) — e(k') 1 is the propaga- 
tor for two non-interacting carriers, and F(k, q, q' ,p, a;) = 
Ei^i^q^Gi^c^/^^S-lFM) is a general- 

ized Green's function related to states where both carri- 
ers have spin up and a magnon is present with momentum 
p. Such states appear when the spin-down carrier flips 
its spin to create a magnon. This process is mediated by 
'H < ex-}x, y and leaves S* ot unchanged. 

We now use Dyson's equation again to obtain the EOM 
for F(k,q,q',p,uj): 



(9) 



lows us to eliminate one of these equations by first 
rewriting Eq. (8) as ^ Eq F(k, q, q' , Q, u)) = -8 qiq > + 
G(k, q, q', ^)G - 1 (t+ 9 , w) - § Eq G(k, Q, g». In- 
sorting this into Eq. (9), taking the sum over p and then 
inserting the result into Eq. (8) we find: 



F(k, q, q',p, u) = { -^L l F ( k > 1, <?', Q, w) - F(k, -q - p, q' , Q, u)\ + J S[G(k, q, g» 



2N 



G(k,~q-p,q',uj)} > G (|+g, f-g-p.W- fl p - J S). 



G(k,q,q',uj) = < 6 q , q > 



1 + y.9(M,^) 



N 



V 



Jo G (±+q,±+Q,il q+ Q) 

2 G (|+Q, %-q,cj- J S) 



2N' 



G 



j w + V — k , k 



G (%+q,%-q,u-J S) 



(10) 



where we used the shorthand Q q = ui — £l q — JqS, and g[k,q,ui) = 1/N Eq Co(|+g, f-+Q, f2 g +Q) can be calcu- 
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latcd numerically. To solve Eq. (10) we use the fact 
that for a chain with N unit cells and periodic boundary 
conditions, the allowed momenta are discrete: g, = j^-i, 
i = 0, . . ., N — 1. Therefore, for any given values of k 
and q' , Eq. (10) is a iV-dimensional system of linear 
equations associated with qi, and we can use basic lin- 
ear algebra algorithms to solve it numerically. Of course, 
one could have used the finite N much sooner and writ- 
ten Eqs. (8) and (9) as two coupled linear systems of 
equations. However, by removing F from Eqs. (8) and 
(9) to obtain Eq. (10), we reduced this to a system with 
N unknowns. Otherwise, we would have to deal with 
an additional N 2 variables for F. This is possible but 
increases the computation time tremendously. 

As already mentioned, model I can be solved similarly 
(for details, see the supplementary material of Ref. 12). 
In both cases, results can be easily obtained for systems 
with N ~ 10 2 , so finite-size effects can be avoided. How- 
ever, note that for these fe-space solutions, it is important 
that N be not too large. The reason is that we are look- 
ing for bound states, and their overlap with the \k,q) 
extended states vanishes like 1/yN, and therefore they 
become "invisible" in the thermodynamic limit. This 
problem is avoided in a real-space formulation, as dis- 
cussed next. 



2. Real-space solution for model I in the two-carrier, 
Siot = NS sector 

This solution is exemplified here for model I. It relics 
on a real-space formulation based on methods introduced 
in Ref. 16. A short discussion of this solution for model 
II is available in the supplementary material of Ref. 12. 

Here we want to calculate the Green's functions: 

G(n,n',k,u) = (k,n'\G(u)\k,n) 
where the states 

ifc(R i+ i+f) 

|jfe,n) = V- f^— ct, 1+C t |FM) (11) 

i 

describe configurations with total momentum k, and 
where the carriers have different spins and are located 
— oo < n < oo sites apart. 

To generate their EOM, we again use Dyson's identity, 
but now we choose % = U +% z s -\-'Hex-,z and V = %—Hq. 
The action of the spin-flip part n£x- x ,y of V links to the 
generalized Green's functions: 

G(n, m, n , k, to) = (A;, n \G(ui)\k, n, m) 

where 



These states describe configurations of total momentum 
k which have both carriers with spin up and located n > 1 
sites apart, plus a magnon located anywhere on its sub- 
lattice, — oo < to < oo. The conservation of S* ot guar- 
antees that states with two or more magnons cannot ap- 
pear, so the EOM involve only these two types of Green's 
functions. These equations are listed in Appendix A. 

To solve this infinite system of coupled equations, we 
use two facts. lh The first is that we can group these states 
and their corresponding Green's functions in terms of an 
integer index M, defined as the distance between the two 
outermost particles. In other words M = n for states 
|fc,n), while for states \k,n,m): 

!n , < to < n 

to , n < to < oo . (12) 

n — to , — oo < to < 

The importance of this index lies in the fact that V 
only connects states for which AM = 0, ±1. We can 
therefore group the 3M + 1 Green's functions with the 
same value of M into a vector V m , and rewrite the equa- 
tions of motions in matrix form: 1 '' 

7mVm = oimVm-x + PaiVm+i, (13) 

where the sparse matrices 7a/, am and /3m can be read 
off directly from the EOM, Eqs. (Al)-(A4). 

To solve this matrix recurrence equation, we use the 
second useful observation, namely that since we are work- 
ing with a finite lifetime I/57, the Green's functions must 
vanish as M — > 00. 1(1 As a result, the solution of Eq. (13) 
must have the form: 

V M =A M V M -i, (14) 

where the continued fraction matrices 

Am = [7m - ftfAif+if 1 a M (15) 

are calculated starting from a cutoff M c for which one 
sets Aj\/ c +i = 0. In practice M c is increased until con- 
vergence within the desired accuracy is reached. While 
this method is used here in ID, it is important to men- 
tion that it works for higher dimensions, as well. 16 As 
an example, in Appendix B we illustrate how to use this 
formalism to calculate G(0,0, k, u>). Other Green's func- 
tions are calculated similarly. 

Of course, we could have used basis states with singlet 
or triplet-like symmetry, such as |fc, n) T ~ \k, n)=F|fc, —n), 
instead of \k,n) of Eq. (11), etc. The usefulness of 
these symmetric states lies in the fact that at k = the 
EOM do not mix the two different symmetries. This re- 
duces the number of equations and therefore computation 
times, and furthermore allows one to identify whether a 
k = eigenstate is singlet or triplet-like. However, for 
any k 7^ 0, the eigenstates do not have a well-defined 
symmetry, i.e. singlet and triplet-like configurations are 
mixed together by the EOM. As a result, if interested 
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in the general solution, there is no advantage in working 
with symmetric basis states. 

For completeness, we also mention that even at k = 0, 
it is incorrect to think of these cigenstatcs as describing 
singlets or triplets of the charge carriers. This is because 
their spins are coupled to the FM background and cannot 
be disentangled from it. Even if the charge-carrier part of 
these states has singlet or triplet-like symmetry, the full 
wavefunction belongs to the Sf ot = NS Hilbert subspace 
and is, therefore, not a proper singlet or triplet. We 
continue to use the terms singlet-like or triplet-like in 
the following purely for convenience. 

Lastly, another advantage of this real-space formula- 
tion is that it allows us to find easily variational solutions 
for the bound states (if any exist), and to understand 
their nature. This is discussed next. 



3. Variational solutions 

The basic idea for the variational solutions (VS) is that 
if bound states exist, their wavefunctions decay exponen- 
tially with the distance n between carriers. We should 
therefore obtain a reasonably accurate approximation for 
them if we neglect all terms in the EOM where the carri- 
ers are farther apart than a certain cutoff; this simplifies 
the calculation tremendously. The VS, however, do not 
predict correctly the continuum states where the carri- 
ers are not bound to each other. This is not a serious 
problem, because the location of such continua can be 
inferred from the one-carrier spectra. 

We considered two examples of VS, with the carriers 
allowed up to 1 and up to 2 sites away from each other, 
respectively. In the following we refer to these as VS1 and 
VS2. As the simplest illustration, we discuss here VS1 
for model I, in the triplet-like sector, and then briefly 
comment on the other cases. 

A triplet-like state allowed in VS1 for model I is: 

e ikRi c 1 . ,|c' i , c: , i . 

i )+ = y V 1 2A h •- u,+ *- t |FM). 

(16) 

Exchange couples it to the one-magnon state: 

p ikRi q— 

|fc,l,l) + = V-^c t 1+ c* 4M f M) (17) 
i v v 

where the magnon is located directly between the two 
carriers. Of course, the magnon is free to move 
so this will further link to the states |fc,l,m) + = 

s, 9iiAit S ^^ 1+m ' FM ) fOT ^ m ^ °- 

Note that the basis state of Eq. (16) is also linked by 
spin- flip processes directly to \k, l,m = 0) + . 

The EOM for the Green's functions associated with 
these (very few) basis states are listed and solved in Ap- 
pendix C. They have 3 poles, i.e. 3 triplet-like bound 



states are predicted by this approximation. Their ener- 
gies (for J — > 0, to simplify their expressions) are: 

w+ = - J J (3 - 4S + Vl + 16S 2 ) , (18) 

u>+ = i J (-3 + AS + Vl + 165 2 ) , (19) 
ujj = 2J S. (20) 

As shown in the next section, the third pole falls inside a 
continuum, in other words it does not survive as a bound 
state if we relax the variational approximation. The other 
two, however, are outside the expected continua and in- 
deed will be shown to give a reasonable approximation for 
the energies of the two triplet-like bound states predicted 
by the exact solutions for this model. 

While this is a nice example of the usefulness of the VS, 
some care is needed. For example, VS1 predicts states 
with triplet-like symmetry for all k. This contradicts 
our previous statement that the exact solution only finds 
states with well defined symmetry at k = 0. The reason 
for this is trivial: no carrier hopping is possible within 
the triplet-like sector of VS1, since there is no "on-site 
triplet", and states with n > 2 arc not allowed. (This 
also explains why the energies of Eqs. (18)-(20) do not 
depend on t or U). This discrepancy is solved by VS2 
where the EOM couple the singlet and triplet-like sectors 
for all k ^ 0, and indeed we find that the bound states 
have well-defined singlct/triplet-likc symmetry only at 
k = 0. This is true for both models I and II. 

VS1 in the singlet sector works similarly, and also pre- 
dicts 3 bound states for model I. Model II is predicted to 
have only one singlet-like bound state and no triplet-like 
bound states. This already suggests that their physics 
may be rather different. Finally, we note that for VS2 
we also restricted the distance between the magnon and 
its closest carrier to be up to 2 sites. This is because 
the magnon propagation is controlled by J, which is the 
smallest energy scale. Further details and more VS re- 
sults arc available in Ref. 18. 



V. RESULTS 

A. Single spin-down charge carrier 

Since the case of a single spin-up carrier is trivial, with 
a single eigenstate of energy E^, here we discuss only 
the single spin-down carrier case. If the spin-flip pro- 
cesses were turned off, %ex\x}y — > 0, eigenstates with a 
total momentum k and S* oi = NS — 5 would either have 
the carrier with spin down and energy Eki, or the carrier 
with spin up and energy Ek—q.\ while the remaining mo- 
mentum and spin is carried by a magnon (spin wave) of 
energy fl q . These latter states give rise to a continuum 
spanning the interval {Ek-qrf + fl q } q . Spin- flip processes 
hybridize the states in the continuum with the discrete 
spin-down carrier state. Of course, the continuum keeps 
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FIG. 3: (Color online) Contour plot of the DOS pi(k, w) for 
(a) model I and (b) model II. The full red lines are the ap- 
proximations given by Eq. (22) and (23) and the dashed 
green lines show the predicted onset of the continuum at 
mm q (Ek- q ,t + fiq). The dotted blue line in panel (a) shows 
the approximation of Eq. (24). Parameters are Jo = 5t, 
J = 0.05i, S = \ and r) = O.Olt in (a), 77 = 0.02i in (b). 



its location, but the discrete state is pushed even lower 
(for Jo > 0) and gives rise to a low-energy, infinitely lived 
quasiparticle: the spin polaron. 

These general expectations are confirmed in Fig. 3, 
which shows contour plots of the density of states (DOS) 
Pl(k 7 uj) = — ^lmGi(k,ui) for models I and II, for J = 5t, 
J = 0.05<, S = |. We only show the lower part of the 
continuum, which indeed begins at mhiq(Ek-q.f + O g ) 
(dashed green line), as expected. The low-energy states 
of both spectra are the discrete states associated with the 
spin-polarons. For Jq < these discrete states appear 
above the continuum, so in that situation the low-energy 
physics is described by incoherent continuum states. 

The nature of the spin-polaron can be understood by 
considering the limit Jo S> t, J. We start with model II. 
The interaction Jiix^ is minimized by an on-site singlet 



between the carrier and its spin: 



1 



V25 + 1 



5" |FM>, (21) 



with eigenenergy — 4^(5 + 1). Strictly speaking, this is a 
singlet only for S = |, but we will use this term in the 

following. Hopping T lifts the degeneracy, and to first 
order in t and J the spin-polaron energy is given by: 



E^\k) 



Jo 



(5 + 1) + 



25 



25 + 1 



e(k) 



2J5 
25+1' 



(22) 



thus 



The spin-polaron state \Pu,k) = -j^J2i e "**" I s , 
describes a singlet, or bound state between the carrier 
and a magnon at the same site, which propagates with an 
effective hopping amplitude 25i/(25+l) < t, suppressed 
because of magnon cloud overlap. The last term is the 
FM exchange energy lost in the magnon's presence. The 
full (red) line in Fig. 3(b) shows this variational approxi- 
mation, which is already quite accurate even for Jq = 5t. 
In fact, Eq. (22) gives a reasonable approximation for 
the k « part of the polaron dispersion down to Jo ~ t, 
showing that this picture of the polaron as a local singlet 
between the carrier and its spin is quite robust. What 
happens as Jo decreases is that the polaron band moves 
closer to the continuum, and the k — > tt part of its dis- 
persion "flattens out" below the continuum, in typical 
polaronic fashion (for more details, see Ref. 14). 

Similar considerations apply to model I. Here, T-L^x is 
minimized by the three-spin polaron (3SP): 19 



\ZSP) 



i+i 



45 



S j ~\~ s, 



i+1 



45+1 



45 



|FM) 



with 1 



It describes a. bound state 



eigenenergy - J (5 + \) 
between the carrier and a magnon on either of the neigh- 
bor spin sites. Hopping lifts the degeneracy and, to first 
order in t and J, the energy of the polaron eigenstatc 

ikR. 



i|35P) 



1+ 



1 is: 



JS 



45+1 45+1 



(23) 



This expression is shown as a full (red) line in Fig. 3(a), 
and is in very good agreement with the exact dispersion 
already for Jo = 5<. The agreement remains reasonable 
down to Jo = t at k w 0, indicating that the description 
of this spin-polaron as a propagating 3SP is also robust. 

For sufficiently large Jo values, model I has a second 
discrete state below the continuum. This is based on an 
excited state of Hex ■ In the Jq ^> t, J limit, it leads to 



\Pf,k) 
energy: 



ikTi 



ASr-Sr )\FM), with 



E^{k) = ^{2S-l)--e{k) + JS. 



(24) 



The agreement between this approximation (thin blue 
line) and the exact result shown in Fig. 3(a) is rather 
poor, especially close to the continuum, but it improves 
with increasing J . More details on these additional 
states are discussed in Ref. 14, in a 2D context. 
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To summarize the results so far, in the Sf, 



NS- 



_ 1 

-'tot IV u 2' 

one-carrier sector, all three models have a low-energy, 
infinitely lived quasiparticle. The mapping from model 
II to model III is straightforward, since the polaron of 
model II has a singlet-like core, precisely like the "hole" 
of model III. Of course, a proper mapping requires some 
rescaling of the parameters: the "hole" hopping integral 
is different from that of the charge carrier itself. The dis- 
cussion above identified the origin of this rescaling (po- 
laron cloud overlap) and the expected renormalization. 

At first sight, the polarons of model I and II seem to 
be quite different. However, we can rewrite 



IP, 



.*> = E- 



JkRi 



f 



4,i, t <sr ) i™) 



i.e. as a singlet between the carrier in an "on-site" orbital 

and its local 



e 2 c: 



2 C 



k,i,<7 V4S+1 V" 

spin. This is exactly like the polaron of model II if we 
replace dk,i,a -> Q t<T . 

(s) 

These dk,i, a operators are similar to the P i o . opera- 
tors in the Zhang-Rice approach: they are linear com- 
binations of carrier orbitals at sites neighboring a spin- 
site, and they are not orthogonal to one another. The 
main difference is that our linear combinations have k- 
dependent phases. This allows us to avoid the need for a 
diverging normalization factor /3k in |P/,fc), which then 
lead to the introduction of the much more complicated 
orbitals <f>i a in the ZRS. 1 While the discussion here is for 
ID, similar considerations hold in 2D. ' ! 

What we found so far, then, is that all three mod- 
els have the same low-energy physics in the non-trivial 
single-carrier sector: the low-energy state is an infinitely 
lived quasiparticle in all cases, with a polaron core 
similar to a ZRS. This suggests that we can indeed 
model the low-energy physics of the two-sublattice, two- 
band Hamiltonian with a much simpler t-J model for 
"holes" , similar to Zhang and Rice's mapping of the two- 
sublattice, three-band Emery model to a t-J model for 
ZRS. The next question is whether this mapping also 
describes well the interactions between quasiparticles. 



B. Two charge carriers in the SI 



NS sector 



As already mentioned, if both carriers have spin up 
(S'tot = NS+ 1), there are no interactions between them. 
If both are injected with spin down (S* ot = NS — 1), the 
magnon-mediated interactions turn out to be weak or 
repulsive and there is no qualitative difference between 
the low-energy spectra of these models, as we discuss at 
the end of this section. 

The interesting case is the S* ot = NS sector, i.e. when 
one carrier is injected with spin up and the other one 
with spin down. As we show now, here we can identify 
qualitative differences between models I and II (model III 
cannot be used here, since it does not distinguish between 
a spin- up carrier and a lattice spin). 
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FIG. 4: (Color online) Spectral function A(k = 0, n = n' = 
l,w) for model I (a) and model II (b). Dashed lines mark the 
expected continuum boundaries. Discrete states, i.e. bipo- 
larons, are labeled according to their symmetry, S for singlet- 
like and T for triplet-like. Parameters are Jo = 20t, J = 0.05i, 
S = 1/2, rj = 0.1* and U/t = 0, 1, 5. 



A(k = 0,n = ri = 1,uj) = -±lmG(n = ri = 
l,k = 0,w) obtained with the real-space solution for 
both models are shown in Fig. 4. Their counterparts 
Aq : q>(k = 0,ui) = — ilmG(fc = 0,q 7 q',uj), are shown in 
Fig. 5. Since both methods are exact, these spectral 
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FIG. 5: (Color online) (a) A q=q i = 2L(k = 0, w) for model I. 
Using q = q' — ^ ensures that states with both symmetries 
have non- vanishing weight, (b) A q = q i =0 (k,ui) for model II 
and k = 0, n, 2n. Dashed lines mark the expected continuum 
boundaries. Discrete states, i.e. bipolarons, are labeled ac- 
cording to their symmetry, S for singlet-like and T for triplet- 
like. Parameters are J /t = 20, J/t = 0.05, S = 1/2, U/t = 0, 
ri/t = 0.02 and N = 60 for (a), N = 50 for (b). 
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weights must have poles at exactly the same eigenener- 
gies, however the weight at these poles differs since we are 
projecting on different types of states. In both cases spec- 
tra are displayed for J$/t = 20, J/t = 0.05 and S = \. 
The value of Jo was chosen so large intentionally in order 
to spread out the different features of the spectrum. 

In the absence of interactions, the spectrum for two 
carriers equals the convolution of the corresponding one- 
particle spectra, which arc known. This is why we ex- 
pect to see a high energy continuum describing scattering 
states consisting of two spin-up carriers and a magnon, 
spanning energies {E q ^ + E q >^ + Q k -q-q'} q ,q' , arising 
from the convolution between the high-energy continuum 
in the spin-down carrier spectrum and the spin-up carrier 
cigenstate. Neglecting a small J term, the edges of this 
continuum arc at ±4t + 2^/JqS, i.e. [16t, 24i] for model 
I and [6t, 14t] for model II, for these parameters. These 
values are marked by the dashed vertical lines. 

The real-space solution indeed shows a high-energy 
continuum for both models at precisely the expected lo- 
cations. However, the k-space solution (Fig. 5(a)) has 
finite weight only in part of the continuum. The rea- 
son for this is that here we project on delocalized states 
where each carrier has a well defined momentum -| ± g. 
Depending on k and q, sections of the continuum have 
vanishing overlap with these states. This is confirmed in 
Fig. 5(b), where A q=q i =0 (k, u) is shown for k = 0,n, 2ir. 
Taken together, the three curves cover the full interval 
where the continuum is expected. Note that the spec- 
trum for k = 2tt is different from that for k = 0, since for 
the former we project on states c\ JFM), whereas for 
the latter we project on cj ^cj JFM). 

The k-space solution also has oscillating weight inside 
the continua. This is a finite-size effect due to the rela- 
tively small value of N = 60 unit cells, and can be fixed 
using a larger N. However, since here we project on delo- 
calized states, an increase in N results in a smaller weight 
for all bound states, as already discussed. Therefore, for 
our purposes it is not beneficial to use a larger N. 

Next, the convolution of the polaron state and the 
spin-up state must result in a low-energy continuum 
{Ep' U \q) + Ek- q ,^} q - Using the asymptotic expres- 
sions from Eqs. (22) and (23) for the polaron energies 
and ignoring the small J term, for k = the edges of 

this continuum are at ±2i 1 1 - 1 - 



0.2 



45+1 



4^ for model I 



and at ±2t ( 1 + 



2S 
2S+1 



ill 

2 



for model II, respectively. 

For the parameters of Figs. 4 and 5 this corresponds to 
[ — 13.7i, — 6.3i] and [— 13i, — 7t], respectively These val- 
ues are marked by vertical dashed lines and are also in 
excellent agreement with the low-energy continua found 
by the exact solutions (similar considerations to those 
discussed above apply to the k-space solution). 

For large values of Jo /t the spin-down carrier spectrum 
of model I has the additional discrete state between the 
polaron band and the continuum (cf. Fig. 3). This 

should also give rise to a continuum {Ep}(q) + Ek- q .f} q , 
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6: (Color online) Exact spectral weight A(k = 0, n = 
n' = l,w) (dashed black line) and the corresponding VSl 
predictions in the triplet-like (top panel, full red line) and 
singlet-like (bottom panel, full blue line), for model I. Param- 
eters are J /t = 20, J/t = 0.05, U/t = 0, S = 1/2, rj/t = 0.1. 



which corresponds to [9t, lit] for these parameters. This 
continuum is absent in the real-space solution. This is not 
surprising because the asymptotic expression for \Pj ,k) 
has zero overlap with any state with a spin-up and a spin- 
down carrier like those used in the projection. In reality, 
there are small corrections to this asymptotic expression 
and indeed, the k-space solution shows a (restricted) con- 
tinuum here, with a weight that is indeed orders of mag- 
nitude smaller than for the other two continua. 

We conclude that there is excellent agreement between 
the continua displayed by the exact solutions, and those 
expected based on the convolutions of the corresponding 
one-carrier spectra; this is a good consistency check. 

If there are sufficiently strong interactions between car- 
riers so as to bind them in pairs, this should result in 
discrete states appearing in the spectrum outside these 
continua. The exact solutions show five such states for 
model I, and one for model II. 

We have checked that these are indeed discrete states 
(as opposed to very narrow continua) by verifying that 
they are Lorentzians of width 77. 18 This finite width is 
responsible for the apparent overlaps of SI and S3 with 
their nearby continua, but they become distinct features 
as 77 — > 0. The nature of these states is further confirmed 
in Fig. 6, where we compare the exact spectral weights 
to the predictions of VSl for the triplet-like (top panel) 
and singlet-like (bottom panel) sectors. Eqs. (18)-(20) 
predict three peaks in the triplet-like sector. The lower 
two are in good agreement with two of the exact peaks; 
the third falls in the center of a continuum, and is there- 
fore not indicative of a true bound-state. In the singlet- 
like sector, VSl predicts three peaks as well. The lowest 
and highest of these happen to fall at the edge of con- 
tinua. A more accurate approximation like VS2 shows 
that they are pushed just outside those continua, where 
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indeed there are discrete states in the exact solution. The 
central peak also agrees well with a discrete peak in the 
exact solution. For model II, VS1 finds one singlet peak, 
in agreement with the peak appearing near lo ~ U in the 
exact solution. 18 

Since the variational solutions include only states 
where the carriers are confined to be close to one an- 
other, the good agreement with the discrete states of 
the exact solutions validates the latters' identification 
as bound states, i.e. bipolarons. Their triplet-like or 
singlet-like symmetry at k = is also in perfect agree- 
ment with what is inferred from the exact solution by 
looking at the relative phases of the n = 1 and n = — 1 
Green's functions, for example. As mentioned, at k ^ 
these states no longer have a definite symmetry. 

From the VS or from the dependence on n of 
G(n, n', k, ui) at energies u> close to one of these peaks, we 
can infer the underlying nature of these polarons. For ex- 
ample, the bound state of model II has most of its weight 
on the configuration where both carriers are on the same 
site, i.e. \k = 0, n = 0). This explains the small spec- 
tral weight for this state in Fig. 4(b), where we project 
on \k = 0, n = 1) states whose contribution is exponen- 
tially smaller, especially for smaller U. It also explains 
the strong dependence of its energy on U : having both 
carriers on the same site costs U. Finally, it also ex- 
plains why the energy of this state hardly depends on Jq: 
the carriers form a singlet with each other when on the 
same site, and there is no exchange between this singlet 
and the lattice spins. As a result, this state of energy 
lu ~ U is always located above the low-energy spin-up 
carricr+polaron continuum, whose energy is ~ — Jq/2 
lower because of exchange between carriers and spins. 
This state (which is very similar to S2 of model I) has 
a non-monotonic dispersion with k in certain regions of 
the parameter space. 18 However, it is rather uninterest- 
ing overall, since it is not bound by magnon-mediated 
interactions and, moreover, it is a high-energy state. 

In contrast, model I has a triplet-like bipolaron (Tl) as 
its low-energy cigenstatc, in addition to four other higher 
energy bipolaron states. While we can analyze and un- 
derstand the nature of all these higher-energy bipolarons, 
in the following we focus on the low-energy bound state 
which is of primary interest to us. 



C. The low-energy bipolaron of model I 

The three configurations with the highest-weight con- 
tributions to this bipolaron's wavefunction are depicted 
in Fig. 7(a). They clearly illustrate the underlying mech- 
anism responsible for binding the two carriers as being 
due to the back-and-forth exchange of a magnon, which 
is emitted by the spin-down carrier and then absorbed 
by the other, initially spin-up carrier. This process is fa- 
cilitated by the fact that the carriers sit on adjacent sites 
and simultaneously interact with the lattice spin located 
between them. The triplet-like symmetry at k = is a 




FIG. 7: (Color online) (a) The three highest-weight configu- 
rations of the low-energy Tl bipolaron of model I, illustrat- 
ing the magnon-mediated interaction that binds the carriers. 
There are matrix elements of order Jo between these states, 
due to carrier-spin spin-flip processes, (b) A possible ana- 
logue in model II for the magnon-exchange process shown in 
(a). The central matrix element is now of order J. Symbols 
have the same meaning as in Fig. 2. 



consequence of the fact that the symmetric combination 
(bonding state) has the lowest energy, all the more so 
since it also avoids the cost of the on-site repulsion U. 
Indeed, Fig. 4(a) confirms that the energy of this state 
is independent of U. 

In Fig. 7(b) we sketch what would be the analog in 
model II for this magnon-mediated interaction. Because 
here each carrier interacts only with its own spin, the 
magnon has to hop from one site to the other. This is 
a process controlled by J, which is the smallest energy 
scale. As such, it is not sufficient to bind the two carriers 
together, and indeed we find no such bound states in the 
low-energy spectrum of model II. Clearly, having carriers 
and spins on different sublattices greatly facilitates the 
magnon-mediated interactions in model I. 

The dispersion E^\{k) of this low-energy bipolaron is 




0.25 0.5 0.75 1 1 2 



k/jc J/t 

FIG. 8: (Color online) (a) The dispersion of the low-energy 
Tl bipolaron of model I, for Jo/t = 15 and J/t = 0.1,0.5, 1. 
(b) Bipolaron mass in units of the free carrier mass, mri/m, 
as a function of J/t for Jo/t — 10, 12.5, 15. Other parame- 
ters are S = |, r)/t = 0.02. Symbols are the exact results 
obtained with the real-space solution, while solid lines give a 
perturbative approximation discussed in the text. 
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FIG. 9: Bipolaron binding energy Eb as a function of Jo for 
two values of J. These exact values are for S = |, and are 
independent of U. The bipolaron unbinds if Jo/t < 6.5t. 



plotted in Fig. 8(a) for three values of J. Changing 
J shifts the curves vertically, as expected since the FM 
exchange energy lost in the presence of the magnon is 
increased. However, the shape of these curves does not 
change much, confirming the fact that the binding and 
dynamics of the bipolaron are primarily controlled by Jo 
and t. The dispersion is consistent with nearest-neighbor 
hopping, but with a renormalizcd effective mass. 

The bipolaron mass expressed in terms of the free car- 
rier mass is shown in Fig. 8(b) vs. J/t, for several val- 
ues of Jo/t. These bipolarons are remarkably light, with 
mri/m ~ 4 for Jo/t = 10 and further decreasing with 
decreasing J and/or increasing J. 

To understand these trends, we used Raylcigh- 
Schrodingcr perturbation theory to third order in t and 
J. 20 Details are summarized in Appendix D. The results 
are shown by the full lines in Fig. 8, and are in reason- 
able agreement with the exact results. In particular, for 
S = i we find the bipolaron effective mass to be: 



171X1 



0.406J 2 /i 
"Jo + 0.570 J 



0.406^ 

t 



0.231-. 
t 



(25) 



The dependence on Jo and J, which agrees with that dis- 
played by the exact results, can be understood as follows. 
Consider moving the 3 configurations shown in Fig. 7(a) 
by one site to the right. For the left and the right config- 
urations, this can be done simply with two consecutive 
hops of the two charge carriers. However, for the cen- 
tral configuration, hopping the first carrier to the right 
costs an energy of order J , since it leaves the vicinity of 
the magnon. Hopping the second carrier then leaves the 
magnon "behind" . This can be fixed either by absorbing 
the magnon first and thus switching to one of the other 
two configurations, or, for large enough J, the magnon 
has a high probability to hop on its own. Of course, in re- 
ality things are more complicated, but this analysis gives 
a rough idea why an increase in J and / or decrease in Jo 
makes the bipolaron more mobile. 

So far, results have been shown for values of Jo/t that 
are rather unphysically large. The final question to ad- 



dress is what happens for smaller values of Jo/t. The 
answer is provided in Fig. 9, where we plot the evolution 
of the bipolaron binding energy Eg (measured with re- 
spect to the lower edge of the low-energy continuum, see 
Figs. 4(a) and 5(a)) with Jo/t. For negative values of 
Eb the bipolaron is below the continuum and therefore 
stable, while Eb — > means that the magnon-mediatcd 
interaction is no longer sufficient to bind the bipolaron 
and instead it dissociates into the continuum of states 
with a polaron plus a spin-up carrier. 

Fig. 9 shows that the bipolaron is stable only for 
Jo/t ^ 6.5, roughly independent of the value of J and 
certainly independent of the value of U. This is consis- 
tent with the fact that Jo controls the magnon-mediatcd 
interaction in model I, and that this interaction must 
be sufficiently large to compensate for the loss of kinetic 
energy in order to bind the carriers. 

An in-depth analysis of the significance of these results 
is offered in the last section. Here we conclude by not- 
ing the extreme importance of the lattice structure to 
the low-energy physics in the two-carrier sector. Because 
in model I the two carriers can simultaneously interact 
with the same spin, the magnon-mediatcd exchange is 
very efficient and results in a strong effective attraction 
that can stabilize a bipolaron. In model II, a spin can 
interact at most with one carrier at a time. As a result, 
the magnon-mediated interaction is very weak and low- 
energy bipolarons do not form. The two models are no 
longer qualitatively similar in this two-carrier sector. 



D. Two spin-down carriers 

It is possible to obtain the exact T = solution in 
the two-carrier, S£ ot = NS — 1 subspace with the real- 
space formalism. The only complication is that spin-flip 
processes link the two spin-down carrier states both to 
states with one magnon and one spin-up carrier, and to 
states with two magnons and both carriers with spin-up. 
Thus, we can now have up to four particles (carriers or 
magnons) in the system, but this can be handled by the 
formalism introduced in Ref. 16. 

However, perturbational arguments described next 
show that it is always energetically unfavorable for a 
low-energy bipolaron to form in this sector. The lowest 
energy feature is, therefore, the two-polaron continuum 
whose location is already known to be {Ep' IT \k — q) + 

E p I,II \q)} q , and the full calculation becomes unneces- 
sary. In model III, the situation is qualitatively similar 
since no interactions appear between two "holes" . 

Consider model I in the limit M0,J^0. If the two 
carriers arc two or more sites apart then each forms a 
3SP, and the total energy is 2E i p I) -> -J (2S+1). If the 
carriers are on neighboring sites, simultaneous interac- 
tions with the central spin change this energy to a value 
that is calculated easily by exact diagonalization. We 
find 18 that this energy is larger than 2Ep^ for any S and 
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Jo, in other words it is energetically favorable to keep 
the 3SP spatially apart. Thus, in model I the magnon- 
mcdiatcd interactions are repulsive in this spin-sector. In 
model II, these interactions are very weak because each 
charge carrier interacts only with its own spin/magnon. 
Turning on the hopping will further favor unbound po- 
larons, since they arc lighter and therefore can lower their 
kinetic energy more than a bipolaron. 

This is why we do not expect low-energy bound states 
in any of these models, in this sector, so we do not con- 
sider it any further. 



VI. SUMMARY AND DISCUSSION 

In this article we used exact, T = solutions for one 
and two carriers doped in a ferromagnetic background to 
compare the low-energy physics of three different models. 
The goal was to understand if the simpler models describe 
the same low-energy physics as the more complex one. 

In the single-carrier sector this turned out to be the 
case: all three models have a low-energy polaron whose 
core can be thought of as a singlet between the carrier 
and a local spin. In the two-sublattice model I, a 3SP 
description is more appropriate but it can be reduced 
to a propagating singlet using a construction somewhat 
similar to the ZRS. Based on this, one could conclude 
that a simple one-band model provides an appropriate 
description for carriers in a FM background. 

However, results from the two-carrier sector paint a dif- 
ferent picture, especially when the carriers are injected 
with opposite spins. Here, for the two-sublattice model 
we find that the ability of both carriers to interact si- 
multaneously with a common spin leads to an attractive, 
magnon- mediated interaction, of order Jo. Such an inter- 
action is evidently not included in a t-J model, where the 
Jo energy scale was integrated out. In fact, for the FM 
background we cannot even use the t-J model to describe 
this case, because "holes" only appear when doping with 
spin-down carriers, and there is no counterpart for a spin- 
up carrier. This is why here we studied the intermediate 
model II, where we know that a spin-down carrier can 
form its singlet-like "hole" , but we can also treat spin-up 
carriers. The fact that in model II only one carrier can 
interact with a given spin, however, leads to a very differ- 
ent mechanism for the magnon-mediated interaction. In 
particular, the strong attraction that appears in model I 
does not arise in model II. 

Similar differences appear in the two spin-down sec- 
tor. We argued that here each carrier forms its polaron 
and they move independently, so the "big picture" seems 
to be the same in all the models: none gives rise to a 
strong-enough magnon-mediated attraction to be able to 
bind the polarons. However, there are differences in de- 
tails. As discussed in the last subsection, for model I here 
we expect a local repulsion between polarons, again con- 
trolled by Jo, due to the additional exchange cost when 
the two spin-down carriers are forced on neighboring sites 




FIG. 10: (Color online) (a), (b) and (c) depict two-hole con- 
figurations on a Cu02 layer, connected to one another by 
spin-flip processes of order Jo. These are the 2D analogs of 
Fig. 7(a), and describe the main mechanism for magnon- 
mediated attraction between holes, (d) and (e) show their 
corresponding one-band configurations, obtained by locking 
the holes in ZRS with a neighbor Cu spin. The thick light 
green lines in (a), (b) and (c) show which hole can pair with 
which spin. For (b), there are two possible pairing options, 
indicated by full and dashed thick green lines, respectively. 



and frustrate each other in forming their 3SP clouds. In 
contrast, in model II each spin-down carrier can form a 
singlet with its spin, and there is no difference if they 
are neighbors or not; hence, there is no counterpart to 
the magnon-mediated repulsion of model I. Model III, 
on the other hand, has only very weak interactions be- 
tween its "holes", due to the fact that if they are on 
neighbor sites they break fewer FM bonds then if they 
are farther apart. This very weak attraction is of order 
J and describes completely different physics. 

Since the magnon-mediated interactions that arise be- 
tween carriers are very different in the three models, it 
is clear that they do not describe the same low-energy 
physics in cases with two or more carriers, even though 
their quasiparticlcs are similar in nature. We therefore 
conclude that the two-sublattice model cannot be re- 
placed by simple one-lattice models, whether with one 
or with more bands, because these no longer allow two 
carriers to interact simultaneously with the same lattice 
spin and, as a result, completely change the nature of the 
magnon-mediated interactions. 

While each new problem has to be considered carefully, 
we expect that this situation may be the rule rather than 
an exception. The reason is that the ZRS-like states on 
which the simpler models are built mix together charge 
and spin degrees of freedom. As such, it seems unlikely to 
us that the magnon-mediated interactions between carri- 
ers could be easily "mimicked" by simple effective inter- 
actions between these complex objects. 

Even though all our exact results are for a FM back- 
ground, we think that simple arguments show that these 
claims extend to AFM backgrounds. Consider, for in- 
stance the situation sketched in Fig. 10(a), with two 
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neighboring holes on a Cu02-like lattice. For simplicity, 
we assume that the Cu spins are well-described by Neel 
order on the length scale of interest here. 

In the two-sublatticc model, configuration (a) is con- 
nected through matrix elements of order Jo to configura- 
tions (b) and (c). This is the 2D analog of the exchange 
of a magnon through the common lattice spin depicted 
in Fig. 7(a), i.e. precisely the process expected to be re- 
sponsible for the magnon-mediated interactions between 
carriers. If we lock these holes into ZRS and map con- 
figuration (a) to a one-band model, we obtain the state 
of Fig. 10(d). Configuration (b) can map to the same 
one-band state (d) , because it is part of the same ZRS as 
(a), but it can also map to configuration (e), depending 
on how we choose which hole pairs with which spins. Fi- 
nally, configuration (c) maps to (e), as expected since it 
is part of the same ZRS as (b). 

The problem is now obvious: the one-band t-J model 
has no matrix element between states (d) and (e), never 
mind one of order Jo (this energy scale has been removed 
from this model) . As such, it clearly fails to describe the 
same physics as the two-sublatticc model, since it grossly 
underestimates the strength of the magnon-mediated at- 
traction between carriers. Instead of this back-and-forth 
exchange of the magnon that strongly favors configura- 
tions with neighbor carriers, the t-J model only describes 
weak interactions of order J between "holes" because of 
the number of AFM bonds broken as they move around 
and reshuffle the lattice spins; this is completely differ- 
ent physics. This example explains why we believe that 
simple one-band models cannot be expected to describe 
properly interactions between carriers. 

Two possible ways to fix the problem can be envi- 
sioned. One is to add additional terms in the one-band 
Hamiltonians to generate the missing matrix elements. 
For example, states (d) and (e) in Fig. 10 could be con- 
nected by a term describing second-nearest neighbor hop- 
ping of a hole if there is a second hole in its vicinity. Sim- 
ilar "conditional" third nearest-neighbor hopping must 
also be included for configurations where the holes are 



on the same line. These terms would strongly favor pair- 
ing of the carriers, especially since they do not disturb 
the AFM background order. In fact, since Jo is a large 
energy scale, this might explain the existence of rather 
mobile pre-formed pairs, which is one of the favorite sce- 
narios to explain the pseudo-gap regime. 

While this way to fix the one-band models is appealing, 
care is needed to make sure that all these additional ma- 
trix elements have the proper signs (i.e., like those arising 
in the more complex model) ; and also that the effects of 
other processes, such as the spin-swap terms'' which we 
have not included in this model, arc properly accounted 
for and "mimicked". It should also be explored what, if 
any, are the consequences of the fact that multiple states 
in the one-band Hilbert space may correspond to a single 
state of the two-sublattice model, as shown in Fig. 10 for 
configurations (b), (d) and (e) (note that there is also a 
third one-band configuration corresponding to (b), with 
the ZRS diagonally opposite to each other). 

The alternative is to study the proper two-sublattice 
models. They have larger Hilbert spaces, unfortunately, 
but at least they are certain to describe correctly the 
significant magnon-mediated attraction between carriers, 
which may well be responsible for most, if not all, the glue 
needed for superconductivity in cuprates. 
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Appendix A: Details for the real-space solution for 
the two-carrier case 

The EOM for model I are: 



1 k 
G(n, n , k,u)) = — < 6 n ,n' — 2i cos — [G(n + 1, n', k, u) + G(n — 1, n', k,u))] 

S 



+ Joy — [G(n, n, n' , fc, tu) + G(n, n + 1, n', k, u>) 



for n > 0, (Al) 



G{n, n , k, u) 



G(0,n',k,ui) 



1 k 

= < 5 n „/ — 2t cos — \G(n + 1, n', k, w) + G(n — 1, n', k, u))] 

lo + ir\ I 2 

s 



~ J o V 2 [G( ~ n ' °' n '' fc ' w) + G( ~"' X ' fc ' w)] | ' 

= \ 5 . n , - 2t cos J [G{l,n', k, u) + G(-l,n', k, u>)} 

uj — U + ii] \ 2 



for n < 0, (A2) 
(A3) 
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G(n, m, n , k, lj) 



lj - 2JS - 4f [45 - (5 m , - 5 m .i - 5 m<n - 5 m ,„+i] + ii] 



x j-t e ''(1 - <5„.i)G(ra - l,m - l,n', fc,w) + e^G(n + l,m + l,n',fc,w) 
+ e _l = G(n + 1, m, n', fc, w) + e l i (1 — 5 n ,i)G(n — l,m,n', k,u) — J5 [G(n, m + 1, n', fc, w) 



+G(n,m- l,n',k,uj)] + J y -[(i m ,„ + <5 m ,„+i)G(n,n', fc.w) 
— (<5m,o + ^m,i)G(— n, n', k, , for n > 1, — oo < to < oo. 



(A4) 



r 



Appendix B: Solution for G(0,0,k,io) 

We show here how to calculate G(0,0, k,uj). From 
Eq. (A3) we see that G(0,0, k,u>) is linked only to 
G(±1,0, k,Lj). We can choose to order the vectors Vm 
in such a way that G(l, 0, k, lj) = [Ai]i,iG(0, 0, k, lj) and 
G(— 1,0, k,oj) = [Ai]2 iG(0, 0, k, lj). Inserting this into 
Eq. (A3) we find: 



G(0,0,k,w) 



1 



in- C/ + 2icosf([Ai]i,i + [Ai] 2 ,i)' 



where Ai is obtained from the recursive relation Eq. 
(15). Once G(0,0, k,Lj) and therefore V is known, 
one can use the recursive relation Eq. (14) to find 
Vm = Am . . . AiV , and thus to generate all Green's 
functions corresponding to values of M > 0. 



Appendix C: Details for VS1 for model I, in the 
triplet-like sector 

We introduce the shorthand notation: 

a+ = + (k,l\G(uj)\k,l) + , (CI) 
a+ A =+<fc,l|GH|fc,l,l) + , (C2) 
at„=+(k,l\G(u)\k,l,m)+ for m < 1, (C3) 



i 1 m =+(fc, l|G(w)|fc, 1, m) + for m < 1, 
and using Eqs (A1)-(A4) we obtain the following EOM: 



1 + J J§ ( a to + V2afi) 

H = " — " , 

lj + irj 

* ia = oj + i v -2JS - J (2S -1)' 

J J§at - JS'V2a+ 1 - J5a+_ x 



l i,o 



w + ir]-2JS- J (25- 1/2) 
J So,^ m j J S 



,m+l 



w + it) - 2 J5 - 2705 



for to < 0. 



(C4) 

(C5) 

(C6) 
(C7) 



As appropriate for VS1, in writing the above equations 
we ignored all Green's functions for basis states where 
the carrier-carrier distance n > 2. 



The last of the above equations is a simple recurrence 
equation that can be solved with the ansatz a^ m (uj) = 
z(ui)a,i m+1 (u)) (note that to < 0). This ansatz makes use 

of the fact that a+ (lj) m ^ °°> and therefore requires 
\z(cv)\ < 1. As previously discussed this is true because 
we are using a finite lifetime 1/rj. 

Then, using af_ 1 (cj) = z(oj)af (uj) in Eq. (C6) re- 
duces the first 3 equations to a linear set with 3 un- 
knowns, 0^,0^ and a^ Q . This can be solved analyti- 
cally. The solution (for J = 0, so as to keep it compact) 
is: 



irj — JqS 



Ld + irj + J (l - 25) 



uj + ir] + J (\ -25) 



(C8) 



While an analytical solution is not possible in all cases, 
the significant reduction in the number of unknowns al- 
lows us to find very easily numerical solutions for the 
various VS, and therefore approximations for the ener- 
gies of the bound states (if any exist). 



Appendix D: Effective mass of the bipolaron 

The first order correction is obtained by diagonalizing 
Hex in the subspace spanned by \k,n = 1) + \k, n = 
— l),\k,n = l,?n = 1) and \k, n = 1, m = 0) + \k, n = 
1, m = 2). Note that all of the above states have triplet- 
like symmetry and that this subspace is invariant un- 
der W^ex since this interaction cannot change the carrier- 
carrier distance. The lowest eigenvalue of T-L^x in this 
subspace gives the zeroth order approximation for Eti, 
while the first order correction in t and J is given by 

the matrix elements of T and Hg for the corresponding 
eigenvector. This results in: 



4°>(fc) 



A4i(Jfc) 



^(3-45+^1 + 1652 



JS 2 + 



5S(-l-4S+v / T+16S 2 ") " 
V1 + 16S 2 

1 + 35 
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Note that (i) the first order correction in t vanishes, 
since T changes the carrier-carrier distance and therefore 
leaves this subspace, and (ii) the zeroth order energy is 
identical to the VS1 result, see Eq. (18), as expected. 

In order to obtain higher order corrections we need to 
include all the states that are linked to by T and l-L^P 
when starting from the configurations considered above. 
If need be, we also have to add whatever other configura- 
tions are necessary to ensure that the enlarged subspace 
is still invariant under T-L^x . This means that we need to 
include a total of seven additional states, some of which 
have singlet-like symmetry, since at finite k mixing be- 
tween the two symmetries occurs. However, for third 



order corrections one does not need to include the state 
where the carriers are on the same site and therefore we 
do not need to worry about U. This may seem counter- 
intuitive, since we are starting out with states where the 
carriers are on adjacent sites, i.e. n = 1, and therefore 
one would expect a transition to n = to be a simple 
process. However, the state with n = only exists for 
singlet-likc symmetry and since we start out with triplet- 
like states, the transition n = 1 — > n = only contributes 
corrections higher than third order. 

For the sake of simplicity we only write here the results 
obtained for S = |: 



a4>) = - ^ 1 '- ■ - 7" 1 — 1 v ~"" 1 (di) 

125(3 + ^) Jo 

(3) (-91 + 125V5) J 3 + 20Ji 2 (-166 + 20V5 + (-209 + 15V5) cos/c) 

T1 W 125CU 2 • (D2) 



(275 + 119V5) J 2 + 40t 2 (245 + 103^ + (205 + 97 V$) cos k) 



Besides an overall shift, a dispersion proportional to cos(fc) arises in the 2nd and 3rd order corrections. The effective 
bipolaron mass can now be extracted easily. 
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